Mean-field regime of trapped dipolar Bose-Einstein condensates in one and two dimensions 
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We derive rigorous one- and two-dimensional mean-field equations for cigar- and pancake- shaped dipolar 
Bose-Einstein condensates with arbitrary polarization angle. We show how the dipolar interaction modifies the 
contact interaction of the strongly confined atoms. In addition, our equations introduce a nonlocal potential, 
which is anisotropic for pancake- shaped condensates. We propose to observe this anisotropy via measurement 
of the condensate aspect ratio. We also derive analytically approximate density profiles from our equations. 
Both the numerical solutions of our reduced mean-field equations and the analytical density profiles agree well 
with numerical solutions of the full Gross-Pitaevskii equation while being more efficient to compute. 



PACS numbers: 03.75.Hh, 75.80.+q, 67.85.-d 

I. INTRODUCTION 

Quantum-degenerate gases with long-range interactions 
have received much attention recently both from experimen- 
tal and theoretical studies. In conventional experiments with 
bosonic quantum gases, short-range interactions have played 
a leading role and are well described by the s-wave scatter- 
ing length [1]. With the realization of a dipolar chromium 
Bose-Einstein condensate (BEC) it is now possible to go be- 
yond such isotropic interactions in degenerate gases [2, 3]. 
Dipolar interactions have a long-range and anisotropic com- 
ponent. These features crucially affect the ground state prop- 
erties [4, 5], stability [6, 7], and dynamics of the gas [8]. Fur- 
thermore, they offer a route for studying exciting many-body 
quantum effects such as a superfluid-crystal quantum phase 
transition [9], supersolids [10] or even topological order [11] 
(for a review of the experimental and theoretical progress in 
dipolar gases see Refs. [12, 13]). 

^^Cr possess a comparatively large magnetic dipole mo- 
ment of 6 Bohr magnetons. A large magnetic dipole moment 
makes atomic BECs an ideal candidate for studying the inter- 
play between contact and dipole-dipole interactions. By re- 
ducing the s-wave scattering length via a Feshbach resonance 
it is even possible that dipole-dipole interactions dominate the 
properties of the BEC [14, 15]. Dipolar effects have also been 
observed in a spinor alkali condensate [16]. Furthermore, 
dc electric fields can induce large electric dipole moments 
in alkali atoms [17]. Systems with a large permanent elec- 
tric dipole moment include heteronuclear molecules [18, 19], 
which are harder to cool to quantum degeneracy [20, 21], and 
Rydberg atoms [22] . 

Complementing the tremendous experimental progress, 
many properties of trapped dipolar BECs have been inves- 
tigated theoretically. The Gross-Pitaevskii equation (GPE) 
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with a nonlocal potential determines its ground state density 
profile in three dimensions and at zero temperature [5]. Ne- 
glecting the kinetic energy term, O'Dell et al [23], Eberlein 
et al. [24] have shown that, remarkably, the ground state den- 
sity profile of a harmonically trapped three-dimensional (3D) 
dipolar BEC is an inverted parabola just as in the nondipo- 
lar BEC. On the other hand, in most experiments to date a 
strong harmonic trap (or optical lattice) along one or two axes 
confines the dipolar BEC to a cigar or pancake shape, respec- 
tively [3, 15]. For these cases, Parker and O'Dell [25] have 
derived one- (ID) and two-dimensional (2D) density profiles 
of a dipolar BEC both in the Thomas-Fermi (TF) limit and 
the ID and 2D mean-field limit. However, their results are 
only valid for polarization along the symmetry axis. Several 
authors have derived effective dipolar potentials in lower di- 
mensions for either axial or transverse polarizations [26-31]. 
To our knowledge, an effective dipolar potential valid for arbi- 
trary polarization direction of quasi- ID and quasi-2D dipolar 
BECs has not been proposed yet. 

Apart from a conceptual clarity, effective equations for 
lower dimensional dipolar BECs also offer a clear advantage 
for numerical computations. For strong trap anisotropics the 
time scales along the compressed and elongated axes are very 
different, which makes an accurate numerical treatment hard. 
Instead of solving the full 3D problem, it is hence desirable to 




FIG. 1. (Color online) In the quasi- ID setup in (a) the dipolar BEC 
is confined to the z direction. In the quasi-2D setup in (b) the atoms 
are confined to the x-y plane. The dipoles are polarized along the 
axis n = {nx,ny,nz) with polar angle (i.e., riz = cos 0). 
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find governing equations for lower dimensional dipolar BECs 
which are suitable for efficient numerical methods. This is 
particularly important for very strong confinement or low den- 
sities, where the usual TF approximation for the full 3D GPE 
becomes invalid [1]. 

In this paper, we present mean-field equations for trapped 
dipolar BECs in one and two dimensions polarized along 
an arbitrary axis. Our equations are based on a mathemati- 
cally rigorous dimension reduction of the 3D Gross-Pitaevskii 
equation (GPE) to lower dimensions. For BECs without dipo- 
lar interactions the formal analysis [32-34] and rigorous anal- 
ysis [35-38] of such a dimension reduction has been discussed 
extensively in the physical and mathematical literature. On 
the other hand, only few mathematically rigorous results are 
available for the dimension reduction of dipolar BECs, such as 
Ref. [39]. The main advantages of our effective equations over 
previously derived results are: (i) valid for arbitrary dipole 
alignment, (ii) well amenable for efficient numerical compu- 
tations typically based on the Fourier transformation and re- 
lated methods. 

For the derivation of the ID (2D) case we assume that the 
BEC is in the ground state of the radial (axial) trap (see Fig. 1). 
We find that the ground state of the lower-dimensional dipolar 
BEC is determined by a modified contact interaction term and 
a nonlocal potential. Compared to a confined BEC with only 
s-wave interactions here the contact interaction also depends 
on the strength of the dipolar interaction and the 2;-component 
of the polarization axis. Crucially, it is independent of the 
transverse components of the polarization. In 2D the nonlocal 
term introduces anisotropy in the BEC ground state density, 
which can be measured, e.g., in time-of-flight experiments as 
a modified aspect ratio of the BEC [3, 40]. We discuss the as- 
pect ratio of the BEC for the intermediate regime between the 
special cases of parallel and orthogonal polarization, which 
have been studied before [23, 25]. Furthermore, we present 
simple analytical density profiles in ID and 2D derived from 
our mean-field equations. We compare ground states of the 
quasi- ID and quasi-2D BEC at different polarizations with the 
ground states of the full 3D BEC and find good agreement. In 
particular, our ground states are a good approximation to the 
ground states of the full GPE in regimes where the TF approx- 
imation fails. In the limit of strong axial confinement we can 
cast our nonlocal potential into a form similar to the Poisson 
equation found for 3D dipolar BECs [23]. 

This paper is organized as follows. In Sec. II, we introduce 
the model of a 3D dipolar BEC at zero temperature. As our 
first main result, in Sec. Ill we present a mean-field equation 
for a quasi- ID dipolar BEC. We compare the ground state so- 
lutions of this ID equation with the full 3D computation and 
an approximate analytical solution. In Sec. IV, we present our 
second main result, namely, a mean-field equation for a quasi- 
2D dipolar BEC. Again we compare its solutions to the 3D 
GPE solution and our analytical approximation. Moreover, 
we calculate the aspect ratio of the BEC if the polarization of 
the dipoles is changed continuously from the longitudinal to a 
radial axis. We conclude in Sec. V. In Appendix A and B we 
present details of the dimension reduction from the 3D GPE 
to the ID and 2D mean-field equations, respectively. In Ap- 



pendix C we derive closed forms for the nonlocal potentials in 
ID and 2D for arbitrary dipole alignment. 



II. 3D MODEL 



We consider a dilute dipolar BEC at zero temperature 
trapped in a harmonic potential V{r) = y (cj^x^ + ^^^^ + 
uj^z'^), where m is the particle mass and ujx,y,z are the trap 
frequencies. We focus on atomic BECs with a magnetic 
dipole moment but it is straightforward to extend the analy- 
sis to degenerate bosonic gases with electric dipole moments. 
We assume that the atoms are polarized along a dipolar axis 



n = {nx.ny.riz) with 
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1. Away from shape reso- 



nances, the wave function V^(r, of the gas is governed by 
the GPE [5, 17, 41] 



'2m 



dd 



V^(r,t), 



(1) 

where g = Airh^as/mi^ the contact interaction strength with 
s-wave scattering length a^. The dipolar potential ^dd is given 
by the convolution 



with the dipole interaction 

Cddl-Scos^O 



Udd{r) = 



47r 



(2) 



(3) 



Here, is the angle between the polarization axis n and the 
relative position of two atoms (i.e., cos 6 = n • r/|r|). For 
magnetic dipoles we have Cdd = MoM^. where /io is the mag- 
netic vacuum permeability and jj^d the dipole moment, and for 
electric dipoles we have Cdd = (P/eo, where eo is the vacuum 
permittivity and d the electric dipole moment. We note that it 
is possible to modify the dipolar interaction Cdd by means of 
a rotating magnetic field [42] . 

We use a mathematical identity to write the dipole interac- 
tion Eq. (3) as [23, 43] 



Uddir) = -Cdd { M ^ 



(4) 



Here we denote with = rixdx^ ^y^y + ^zdz the derivative 
along the dipole axis and dnn = Sn(^n)- Inserting Eq. (4) and 
Eq. (2) into the GPE, Eq. (1), results in 



2m 



(5) 

We note that two terms contribute to the dipolar interac- 
tion. The first term in Eq. (4) reduces the contact in- 
teraction strength [third term in Eq. (5)], while the sec- 
ond term in Eq. (4) contributes the potential = 

-CrfrfSnn /d3r'/73D(r-rO I V^(r^t)p with kernel /73D(r) = 
l/47r|r|. 
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We introduce dimensionless quantities by rescaling lengths 
with r rao, times with t t/uo, energies with huo, and 
the wave function with ifj^N/aQ, where ujq is the small- 
est trap frequency in the system [ujq = min(cja^, cj^, cj^)], 
ao = yJh/mujQ the corresponding magnetic length, and N 
the total number of atoms in the BEC. After rescaling, in di- 
mensionless form Eq. (5) is given by 



-V^ + F(r) + /^(l-edd)|V^|' + ^3D 



V^, (6a) 



^3D = -SpCdddnnJ d^r'U3j,{r - r')\^{r\t)\^ (6b) 

where P = AnNas/ao and Cdd = Cdd/^9 defines a natural 
dimensionless parameter for the relative strength of dipolar 
and s-wave interactions. The dimensionless trapping potential 

is V{r) = \{nlx^ + -ilv^ + -iiz^) with 7^ = cj^/cjq, ly = 



III. QUASI-ID DIPOLAR BEC 

By choosing a sufficiently large radial trap frequency it is 
possible to freeze the radial motion of the BEC [44]. If the 
extent of the radial cloud is much larger than the s-wave scat- 
tering length, this is the limit of a quasi- ID BEC [45]. In 
Fig. 1(a) we illustrate the geometry of this setup. In this sec- 



tion, we present an intuitive mean-field equation for the axial 
wave function of such a strongly confined dipolar BEC. Our 
equation is based on a reduction of the 3D GPE to ID assum- 
ing a strong radial confinement. 



A. ID mean-field equation 

In order to derive a mean-field equation for the axial wave 
function of the condensate we assume that Uz <^ ujx = ujy =: 
uj± and gno <C huj_\_, where no is the peak density of the BEC. 
Moreover, we require that dipole interactions do not excite 
radial modes, i.e., ^dd ^ huj±. Then the radial modes of the 
BEC are in the ground state of the transversal harmonic trap 
and the order parameter of the BEC factorizes. We write 
the factorized wave function as 

^{r,t) = e-'^^'w2D{x,y)^ij,{z,t), (7) 
W2Mx^y) = ^l^e—^^'^y'y'\ (8) 

In this section we rescale equations in terms of the dimen- 
sionless lengths r ra^, times t ^ t/ujz, and the axial wave 
function ?/^id i^iD\/ N/az with = ^JhJmhTz the mag- 
netic length in z direction. Energies are expressed in units of 
hijJz. Given these assumptions, in Appendix A we show that 
the 3D GPE, Eq. (1), reduces to an equation for the axial wave 
function V^id 



idtipiD{z,t) 

^ID 



{-^S.. + ^id(^) + ^ fl + ^ (1 - 3n^)l |V^id(^, t)p + \ M^.t), 



(1 - 3n^) dz. r dz'U,^{z - z')\M^'M 

oyZTT J-oo 



(9a) 
(9b) 



where Vixy{z) = z'^/2, /3id = gN/huJzCtl = 4:7rNas/az and 
riz is the z component of the dipole axis n. The trap aspect 
ratio is given by j = uj±/uJz- We find for the kernel /7id 

Ui^{z)=e'''"/^evk{\z\./^), (10) 

where erfc is the complementary error function. In Ap- 
pendix C we calculate the derivative in Eq. (9b) and give a 
closed form of the resulting convolution integral. However 
for our focus on numerical computation of the ground state, 
expression (9b) is better suited. It is worth pointing out that 
the only approximation in the derivation of Eq. (9) from the 
GPE is the factorization in Eq. (7) with the choice of a Gaus- 
sian as the radial wave function W2d{x, y). 

The formula for the wave function of a quasi- ID dipolar 
BEC, Eq. (9a), is very intuitive. It has the same structure as 
the corresponding 3D expression Eq. (6a). Notably, the effect 
of the dipolar interaction in ID is an altered contact interac- 
tion strength and the introduction of a nonlocal potential <I>id. 
In Fig. 2(a) we plot the kernel /7id of the potential, Eq. (10), 



for different trap aspect ratios. In contrast to the kernel [/sd in 
3D, the kernel of the ID potential does not diverge at the ori- 
gin. Instead we find /7id(^) = 1 - v^kl/v^ + 0{\z\^) 
for |z| 0. In the opposite limit, \z\ oo, we find 
^iD (^) ~ >/2/ I ' ^- ' it scales with l/\z\ as the 3D ker- 
nel. These properties allow for efficient numerical methods 
based on Eq. (9) [43]. 



B. Ground state 

The third term in Eq. (9a) describes the altered contact 
interaction, which now depends on the dipolar interaction 
strength Pw^dd- The anisotropy of the dipolar interaction is 
manifest in this term in the dependence on the z component 
Uz = cos of the dipole axis, where is the angle between 
dipole and z axis [cf. Fig. 1(a)]. If the dipoles are aligned 
along the longitudinal BEC axis (Uz = 1), the effective con- 
tact interaction reduces by a factor of{l — edd)- Neglecting the 
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FIG. 2. (a) The kernel Um, Eq. (10), for 7 = 10 (solid line), 80 
(dashed), 1000 (dotted), (b) Uid, Eq. (17), for 7 = 1/10 (solid), 
1/80 (dashed), 1/1000 (dotted). In both cases increasing the con- 
finement of the BEC leads to increasingly local behavior. 



nonlocal part of Eq. (9a) this results in a reduced BEC length. 
If the dipole axis is perpendicular to the BEC axis (n^ = 0), 
the BEC length increases since the contact interaction is larger 
by a factor (1 + 6^^(^/2). Intuitively, these two cases can be un- 
derstood in terms of a string of magnets: magnets with poles 
aligned head to tail attract each other, while magnets in a head 
to head or tail to tail configuration repel each other. However, 
our equation also maps out all intermediate configurations be- 
tween these two special cases. Finally, we notice the familiar 
increase of the effective s-wave scattering strength in ID by a 
factor 7/27r, which is due to the choice of a Gaussian ground 
state in the factorization Eq. (7) [33, 46]. 

The potential <I>id in Eq. (9a) describes the nonlocal effect 
of the dipolar interaction on the BEC. The shape of the kernel 
[/id [see Fig. 2(a)] reveals that this potential becomes more 
local with increasing trap aspect ratio 7. Moreover, for large 
7 we expect that <I>id does not affect the shape of the BEC 
significantly. Owing to the properties of convolutions we may 
apply the second derivative in <I>id only to the density |^idP 
in the integral. However, for the ground state we expect this 
density to become flatter for large 7 so that the derivative be- 
comes smaller. Since the contact interaction scales linearly 
with 7, the contact term dominates over the nonlocal poten- 
tial. Similar to the modified contact interaction our Eq. (9b) 
explicitly states the dependence of <I>id on the z-component 
of the dipole axis (and predicts no dependence on other com- 
ponents). 

We have seen that the dipolar interaction in ID is com- 
posed of a local (or contact) interaction and a nonlo- 
cal interaction. In order to determine the sign of the 



nonlocal interaction, we evaluate the corresponding en- 
ergy via Fourier transformation. The energy is given 

by ^fdz^ij)\iljiB{z,t)\^ = ^ /c^/c^S^IV^iDp {k^,t) = 



3/3iDeddV^ /-i 
16 V-"-" 



, where 



denotes the Fourier trans- 



f{h) = {l/V2^)Jdzf{z)e- _ 
form of a function /. Since the Fourier transforms /7id and 
I IV^iDp I are positive (see Appendix A), the sign of the non- 
local potential is opposite the sign of the modification in the 
contact term ^^^^(1 — 3n^). This means that the nonlo- 
cal potential counteracts the action of the contact term: if the 
contact term is repulsive, the nonlocal potential becomes at- 
tractive and vice versa. 

Another useful observation is the following. The BEC 
ground state is not affected by the j)resence of dipolar inter- 
actions for the "magic angle" cos 6m = l/V^ [47]. In this 
case, the 3D dipolar potential Udd vanishes and accordingly 
Eq. (9) reduces to the GPE for a quasi- ID BEC without dipo- 
lar interaction. This observation could be useful for very sen- 
sitive matter wave interferometers, where the dipole interac- 
tion dominates the decoherence when the s-wave scattering 
length has been reduced via a Feshbach resonance [15]. 

We can derive an analytical solution for the density of 
the quasi- ID dipolar BEC if the contact interaction term in 
Eq. (9a) is repulsive and dominates the dynamics [1]. In this 
case, we neglect the kinetic and nonlocal parts in Eq. (9a). 
Assuming a stationary solution, V^id(^, t) = e~*^^^ Y^niD(^), 
we find the density profile niD(^) = [/i^ — (2;^/2)]{/3id 2^ [1 + 



^(1 



3n^)]}~^, where /i^ is the chemical potential along 
z. The density vanishes for z > Z = \^2jj^. We obtain the 
halfwidth of the condensate Z by evaluating the normalization 
condition J dzn{z) = 1, which yields 



'3Ad7 
2 27r 



l + ^(l-3n^) 



1/3 



Inserting this expression into the density profile gives 

3 / .2 



nw{z) 



4Z 



1 



Z' 
Z2 



(11) 



(12) 



Parker and O'Dell [25] have derived analytical ID densities 
for the special case of dipoles aligned along the z axis. In 
this case, we note that our generalized expression Eq. (12) 
coincides with their density in the limit termed "ID mean- 
field regime". Generalizing the criterion for the validity of the 
mean-field regime in Refs. [25, 48], we see that Eq. (12) is a 
good approximation for the axial density profile of an elon- 
gated dipolar BEC if 
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ID 



^dd 



l + ^(l-3n^) 



< 1. 



(13) 



In Fig. 3 we compare the density of the quasi- ID BEC de- 
termined via Eq. (9) with the analytical prediction Eq. (12) 
and the numerical solution for the full 3D GPE in Eq. (1) af- 
ter integrating over the transversal directions. If the dipole 
axis points perpendicular to the elongated BEC axis, we can- 
not distinguish the analytical density from the the quasi- ID 
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FIG. 3. (Color online) Linear density of the quasi- ID BEC according 
to the solution of our 2D equation, Eq. (9) (blue solid lines), the 
corresponding analytical prediction of Eq. (12) (gray dotted), and the 
full 3D GPE of Eq. (1) (shaded area). In the upper panel dipoles are 
aligned with the BEC axis, while in the lower panel they are aligned 
perpendicular to the BEC axis. We choose /^m = 100, tdd — 0.9 
and the 7 given in the plots. 

GPE result (see lower panel of Fig. 3). We have checked that 
for larger values of the trap aspect ratio 7 the discrepancy be- 
tween the full GPE result and our approximate ID solution 
with perpendicular dipoles diminishes. On the other hand, if 
the dipoles are aligned with the BEC axis our ID mean-field 
solution agrees very well with the solution of the full GPE (see 



upper panel of Fig. 3). Furthermore, we note that the BEC is 
compressed compared to the case with perpendicular dipole 
axis. This compression is caused by the attraction of aligned 
dipoles in a ID setup, and it is manifest in the reduced con- 
tact interaction in Eq. (9a) for = 1. However, for axial 
dipoles the nonlocal term in Eq. (9a) produces an appreciable 
repulsive potential. As a consequence, the BEC is broadened 
as compared to the analytical profile. Therefore, in the regime 
of small or moderate interaction energy /3id, we see that the 
usual approach of neglecting the kinetic and nonlocal terms 
is not sufficient to describe the density profile. On the other 
hand, our proposed ID equation, Eq. (9), describes the BEC 
accurately in the mean-field regime at experimentally relevant 
trap aspect ratios 7. 



IV. QUASI-2D DIPOLAR BEC 

In this section we consider a dipolar BEC which is strongly 
confined along the z axis, 7 <C 1 [cf. Fig. 1(b)]. Anal- 
ogous to the preceding section, we assume that gn^ <C huj^ 
and that the axial extend of the cloud is much larger than the 
s-wave scattering length. If the dipolar interactions are also 
small compared to the axial trap energy hjj^, then the BEC is 
in the ground state of the axial harmonic trap. This is the case 
of a quasi-2D BEC, where the BEC wave function separates 
into 

V^(r, t) = e-'^^'^^^2D{x, y, t)wM, (14) 

^^^(,)=f!!^)^/\--^.^V2^ (15) 
v Tih J 

In this section we use the dimensionless rescaling r ra^, 
t t/uj^, ?/^2D i^TD \/N/a\, where a±_ = ^/h/rruJ±, is 
the radial magnetic length. Energies are given in units of huj±_ . 
We show in Appendix B that the transversal wave function 
ijj2D fulfills the following equation 



^l)TD{x,y,t), (16a) 



2D 



[Sn^n^-n^V'j j dx'dy'U2D{x-x\y-y')\^2D{x\y\t)\^' 



(16b) 



Here, V2D{x,y) = (x^ + y'^)/2 and (32d = ^7rNas/a±. We 
denote with n^dx + riydy and dn^n^ dn^{dn^). 

The kernel U2D is radially symmetric and is given by 

where Kj^ (u real) denotes a modified Bessel function of the 
second kind and = {x — x')'^ ^ {y — y')'^ . In Appendix C 
we show that the nonlocal potential, Eq. (16b), can be written 
as a simple convolution <I>2d = — J dx' dy' [x — 



x' ^y — y')\ iIj2d {x' ^y' ^t)\^. There we also derive a closed form 
for , which explicitly depends on the polarization axis. 
Assuming validity of the GPE, the only approximation in the 
derivation of Eq. (16) is the factorization Eq. (14). 

In Fig. 2(b) we plot the kernel /72d, Eq. (17), for different 
trap anisotropics 7. In contrast to the equivalent ID kernel 
/7iD in Fig. 2(a), the long range behavior does not depend on 
7. In fact, we can show that /72d(^) ^ l/27rr for r 00. 
This is equivalent to the long-range behavior of the 3D kernel 
/73D. However, in the opposite limit, r ^ 0, we find that 
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the divergence of the kernel is only logarithmic, ^72d(^) 
ln(r) + ln(2v^) + const]. 



For numerical computations in Fourier space, the expres- 
sion Eq. (16b) for the nonlocal potential is often more useful 
than the closed form derived in Appendix C. Moreover, in the 
limit of large trap anisotropy, 7 <C 1 we have shown that the 
potential in Eq. (16b) is equivalent to a Poisson-type equation. 
To this end we introduce the fictitious potential ^2D defined by 
^2D = -^^^[Sn^n^-n^V2]^2D. For 7 < 1 wc may then 
replace Eq. (16b) by 



(18) 



Hence, the computation of the nonlocal potential ^2d in 
Fourier space involves only multiplications of the density 
|V^2Dp with the momentum. 

In contrast to the ID mean-field equation, Eq. (9), in 2D the 
dipolar interaction increases the contact interaction strength 
for dipoles aligned along the z axis (and positive Cdd)- This 
is a manifestation of the fact that magnets aligned in parallel 
repel each other. The modification of the contact interaction 
term by a factor of 1/ ^/2iry is due to the compression along 
the z axis [33, 45]. Furthermore, unlike in ID the effect of 
dipolar interactions does not vanish at the magic angle Om'- 
while the dipolar contact interaction term vanishes, the nonlo- 
cal term [last term in Eq. (16a)] does not. 

Analogous to the quasi- ID case, we now derive an approx- 
imate analytical expression for the density. To this end we 
assume that a repulsive contact interaction term (third term) 
in Eq. (16a) dominates the ground state solution. Hence, we 
neglect the kinetic and nonlocal terms in Eq. (16a). With the 

-2/xrty/^2D(r) we find the 



stationary ansatz ^2d(^, ^, t) 



e 

density profile nTD{r) = [fir - {r^ /2)]{I32Di- ^2¥^v 
3n^)]}~^, where = + y'^ and /i^ is the radial part 
of the chemical potential. The density vanishes for r > 
R = \J2[ij.. By evaluating the normalization of the density, 
27r / drrn2D{r) = 1, we find the mean-field radius 



[1 



;(1 



R = 



4/3: 



'2D 



[1 - edd(l - 3n^ 



1/4 



(19) 



Inserting this radius into the analytical density profile yields 



n2D{r) 



ttR^ 



1 



(20) 



For the special case = 1 our expression for the density 
^2d(^) coincides with the expression for the "2D mean-field 
regime" given in Ref. [25]. We can formally generalize the 
condition for the validity of Eq. (20) given in Refs. [25, 48] to 



^2DV7 

47r 



(21) 



While this condition may suggest that n2D(^) is a good ap- 
proximation for large dipole moment and small axial polariza- 
tion Uz, we note that in the regime ^ 1/ a/3 the anisotropy 
and magnitude of the potential ^2D increases appreciable. 



This may be seen by evaluating the kernel given in Ap- 
pendix C. In other words, n2D{r) is a good approximation for 
the 2D density profile if condition Eq. (21) is fulfilled and the 
dipoles are polarized predominantly in the axial direction. 

By numerically solving Eq. (16) we obtain radial den- 
sity profiles of a quasi-2D dipolar BEC for various trap 
anisotropics and polarizations. For axial polarization (n^ = 1) 
we find a radially symmetric density. For nonaxial polariza- 
tion (riz < 1) we find that the quasi-2D BEC is elongated 
along the polarization axis projected onto the x-y plane and 
compressed orthogonal to the polarization axis. This is in con- 
trast to the quasi- ID case where the attraction between aligned 
dipoles compresses the BEC along the polarization axis. This 
is a result of the saddle shape of the potential with min- 
ima along the projection of the dipole axis. The experiments 
in Ref. [3, 49] show such an elongation of a dipolar BEC in 
the Thomas-Fermi regime. 

In Fig. 4 we show density profiles of the quasi-2D BEC 
along the elongated (solid blue lines) and compressed axes 
(dashed red) in the x-y plane. If the dipoles are aligned par- 
allel to the symmetry axis of the quasi-2D BEC (n^ = 1), 
the overall dipolar interaction is repulsive. This is manifest in 
Eq. (16), where the contact interaction term becomes larger 
for positive dipole strength Cdd and the nonlocal potential is 
positive. Moreover, the BEC remains radially symmetric as a 
result of the vanishing radially asymmetric derivative dn^n± • 
We plot the radially symmetric density profile for = 1 in 
the top panel of Fig. 4. The radial profile of the BEC be- 
comes increasingly asymmetric as we move the polarization 
away from the z axis. This is evident in the different widths 
of the density profiles along the two orthogonal axes in the 
middle and bottom panels of Fig. 4. The plots with polariza- 
tions along the x axis or the diagonal of the x-y plane (mid- 
dle and bottom left panel of Fig. 4) show the largest differ- 
ence in width. The case of equal polarization in all directions 
(bottom right panel in Fig. 4), n = ^(1, 1, 1), is special 
because the dipole interaction does not have a local charac- 
ter. This is manifest in Eq. (16) where the contribution of 
the dipole interaction to the contact interaction vanishes at the 
angle = cosO ^ 54.7°. Consequently, in the bottom right 
panel of Fig. 4 we only observe a very small asymmetry of the 
radial BEC density, which is a purely nonlocal effect caused 
by the potential ^2D- 

For comparison, in Fig. 4 we also plot the density pro- 
files of a dipolar BEC obtained by numerically solving the 
3D GPE, Eq. (1), and integrating over the z direction. We 
find excellent agreement with the solutions of our proposed 
2D equations (16) for sufficiently large trap anisotropics. In 
the top left panel of Fig. 4 we notice a slight discrepancy to 
the 2D solution because the trap anisotropy is not sufficient 
to suppress the change in the axial density profile caused by 
the dipole interaction. Furthermore, in Fig. 4 we plot the ap- 
proximation n2D{r), Eq. (20). We observe good agreement 
within its regime of validity discussed below Eq. (21). The 
analytical approximation of the radial profile agrees well with 
the numerical ground state solution if the polarization is pre- 
dominantly perpendicular to the BEC disc (top and bottom 
right panels of Fig. 4). For polarizations predominantly in the 
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FIG. 4. (Color online) Cuts through the radial density profiles of the 
quasi-2D dipolar BEC given by Eq. (16) for various polarizations and 
trap anisotropies. The cuts are taken along the axes with largest (x 
axis, solid blue lines) and smallest extend of the BEC (y axis, dashed 
red). The insets show density plots of the quasi-2D BEC and the lines 
indicate the position of the cuts (x and y axes, respectively). The gray 
dotted lines are the analytical profiles n2D(r) and the shaded areas 
are the profiles obtained from the 3D GPE, Eq. (1). For sufficiently 
large confinement the 3D GPE profiles are not distinguishable from 
our 2D solution. We choose (Sid = 100, Cdd = 0.9 and the dipole 
axis n = (0,0, 1) (top panel), n = (1,0,0) (middle panel), n = 
^(1, 1, 0) (bottom left panel) and n = ^(1, 1, 1) (bottom right 
panel). 



plane of the quasi-2D BEC our proposed Eq. (16) remains a 
good approximation while the analytical approximation be- 
comes invalid. 

We have seen in Fig. 4 that the quasi-2D BEC in a radi- 
ally symmetric trap loses its radial symmetry if the dipole axis 
does not point along its symmetry axis. This is a consequence 
of the anisotropic nature of the dipolar interaction. It is pos- 
sible to observe this effect experimentally by measuring the 
aspect ratio 



ay 
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FIG. 5. (Color online) Aspect ratio of the quasi-2D BEC for (a) 
varying dipole strength Cdd with polarization along the x axis and (b) 
varying polarization angle in the x-z plane [n = (sin^, 0, cos^)] 
with Cdd = 0.9. We use the trap aspect ratios 7 = 1/10 (solid 
lines) and 7 = 1/80 (dotted) with /32d = 100. The circles indicate 
the corresponding condensate aspect ratio according to the numerical 
solution the 3D GPE, Eq. (1). The upper axis in (b) shows Uz = 
cos 0. 



of the BEC. The aspect ratio is particularly suited for mea- 
surements because it is not sensitive on the exact number of 
particles [40]. 

We have computed the aspect ratio of a dipolar BEC by nu- 
merically solving our 2D equations (16) for different values 
of the dipole moment and polarization angle. Fig. 5 summa- 
rizes these results. In Fig. 5(a) we plot the aspect ratio for 
polarization along the x axis and varying dipole interaction 
strength Cdd- We observe that the BEC becomes increasingly 
elongated as we increase the dipolar interaction strength. In 
Fig. 5(b) we fix the dipole strength but let the dipolar axis ro- 
tate in the x-z plane. This shows how the radial BEC density 
profile changes from a symmetric disc to an oval shape as we 
rotate the polarization away from the z axis. We note that 
large changes in the aspect ratio occur in a region around the 
magic angle 6 6m, where the effective contact interaction 
is nearly independent of the dipole interaction strength Cdd- 
Large trap anisotropies [dashed line in Fig. 5(b)] suppress the 
onset of significant BEC asymmetry because the contact in- 
teraction dominates the ground state. We have also obtained 
the aspect ratios by numerically solving the 3D GPE, Eq. (1) 
and integrating over the z axis (circles in Fig. 5). We find 
excellent agreement with the aspect ratios obtained from our 
2D equations. Since a rotation of the polarization in the x- 
y plane only corresponds to a rotation of the elongated axis 
of the ground state density (see Fig. 4), we obtain similar re- 
sults to Fig. 5(b) for arbitrary polarization. This shows that 
the reduced Eqs. (16) are indeed a good approximation for de- 
scribing quasi-2D dipolar BECs at arbitrary polarization with 
sufficiently strong axial trapping. 



V. CONCLUSION 

We have presented Gross-Pitaevskii-type mean-field equa- 
tions for trapped quasi- ID, Eq. (9), and quasi-2D, Eq. (16), 
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dipolar BECs polarized along an arbitrary axis. These equa- 
tions are based on a rigorous dimension reduction of the full 
3D GPE. In contrast to previous works, they are valid for ar- 
bitrary dipole alignment in the mean-field regime if the BEC 
is in the ground state of the radial or axial harmonic trap, re- 
spectively. Our result shows that quasi- ID and quasi-2D dipo- 
lar BECs are governed by a modified contact interaction term 
and an additional nonlocal potential. We have given explicit 
expressions for the nonlocal potential for arbitrary polariza- 
tion (also see Appendix C). 

One of the main advantages of the proposed mean-field 
equations is that they are well suited for numerical compu- 
tations in strongly confined BECs. Our numerical implemen- 
tations of the ground state computation in ID and 2D perform 
much faster than our equivalent 3D computations. This is be- 
cause we only need to integrate over the reduced dimensions, 
which vary over similar time scales, whereas the excluded 
dimensions vary on a much faster time scale for strong trap 
anisotropics. Moreover, the kernel of the convolution in the 
nonlocal potential is bounded in ID and diverges only loga- 
rithmically in 2D. In contrast, in 3D the corresponding kernel 
diverges inverse linearly. Our formulation of the nonlocal po- 
tential in terms of partial derivatives allows for efficient nu- 
merical methods based on the Fourier transformation. 

We have computed the ground states of our ID and 2D 
equations numerically and compared them with the ground 
states of the 3D GPE. We find excellent agreement but no- 
tice small discrepancies in ID for the case when the BEC 
is polarized perpendicular to the elongated direction. By 
neglecting the kinetic energy term and assuming a vanish- 
ing nonlocal potential we have derived analytical expressions 
for the density profiles of quasi- ID and quasi-2D dipolar 
BECs when the dipoles are aligned predominantly along the z 
axis. The ground state of the quasi-2D dipolar BEC becomes 
anisotropic if the polarization is not parallel to the symmetry 
axis of the pancake-shaped BEC. This results in a varying as- 
pect ratio of the BEC pancake, which can be measured with 
time-of-flight imaging. We have computed this aspect ratio for 
varying dipole interaction strength, trap anisotropy and polar- 
ization axis. 

As an outlook, we point out that it is straightforward to use 
our equations for studying the dynamics of arbitrarily polar- 
ized dipolar BECs in lower dimensions. They also allow for a 
stability analysis of such lower dimensional dipolar BECs. In 
a future work we will investigate the influence of a changing 
dipole axis on the formation of vortices in 2D dipolar BECs 
in a rotating frame. 
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Appendix A: Derivation of the ID mean-field equation 



Under the assumption in Sec. Ill, plugging V^(r, t) = 
e~'^^^^W2D{x, y)ilJiB{z^ t) into the GPE, Eq. (6a), after rescal- 
ing, we have 



(Al) 



where 



^ = -^x^edd j d^rU3D{r - rO|(V^iD^2D)(r , t)!^ (A2) 



with the rescaled quantities W2Y){x^y) = y^e~i^^ +^ \ 
Pi = /^id(1 — ^dd)- By multiplying both sides of Eq. (Al) 
by W2d{x^ y) and integrating over the x-y plane we obtain an 
equation in z only. Since W2t>{x^ y) is normalized according 
to J dxdyw2j){x^y) = 1 and J dxdyw\jy{x^y) = 7/27r, this 
integration is straightforward for all but the last term. In the 
following we outline the calculation of this last term. 
In order to obtain Eq. (9), we need to compute 



dxdydnn^{x, y, z, t)w2j){x, y). (A3) 



Noticing the symmetry of U^d and W2d in x and y we see that 

dxx {U3D * |V^lD^2Dp) = dyy {U3D * |V^lD^2Dp) , (A4) 

where we denote with a * 6 the convolution of a and b. Since 
U^j) is the Green's function of the Poisson equation [23], we 
have — V^/73D(r) = S{t). Inserting Eq. (A4) into this expres- 
sion yields 



dxx {UsD * |^1D^2d|^) 



|V^lD^2Dp 

2 

dzz {UsD * |V^1D^2d| 



(A5) 



Moreover, |V^id^2dP is an even function with respect to x 
and y and so is U^j) * | V^id^2D which implies that the partial 
derivatives dxy, dyz, dxz of Usj) * |V^id^2dP are odd functions 
in X and y. Recalling d^n = nld^x + '^^l^yy + nld^z + 
2nxnydxy + 2nynzdyz + 2nxnzdxz we plug Eq. (A5) into 
Eq. (A3), which yields 



27r' 



3ni 



J J dxdy {U3D * |V^iD^2Dp) ^Id^ (A6) 



because integrals of odd functions (the mixed partial deriva- 
tives) vanish. We see that the first term in Eq. (A6) contributes 
to the modification of the contact interaction term in Eq. (9a). 
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In order to evaluate the last term in Eq. (A6) we expand the 
convolution and compute the resulting integrals 



R4 47Ty^{x - X'Y + - y'Y + (2: - Z'Y 



To this end, we introduce the variables x = x—x' ,x' = x-\-x\ 
y = y — y' ,y' = y^y' and integrate over x' and y' . The result 
is 



dxdy 



76 



87r2 ypTFTF^^^ 



(A7) 



which can be evaluated further in polar coordinates x -- 
r cos 0',y = r sin 0'. Equation (A7) reduces to ^^^ ?7id(^ 
z') with the substitution ( = y/r'^ -\- {z — z'Y and 



TT 



dC,e 



-7CV2 



(A8) 



\A 



Equation (A8) is an integral representation of /7id given in 
Eq. (10). Inserting Eq. (A8) back into Eq. (A6) we obtain 
the potential $id, Eq. (9b). Plugging Eq. (A6) into the inte- 
grated Eq. (Al) results in the mean-field equation for a quasi- 
ID dipolar BEC, Eq. (9). 
The Fourier transform of t/iD(^) is given by 



f/lD(^.) 



1 



e 27 



(A9) 



The asymptotic behavior of this Fourier transform is 

U^i^z) = ^[ln(27)-7e-21n|fc,|] + 0(A:,)for|A:,| ^0 

and /7id(^2) ^ l^^l ^ where 7e is the 

Euler-Mascheroni constant. 



Appendix B: Derivation of the 2D mean-field equation 

Under the assumptions in Sec. IV, plugging V^(r, t) = 
g-iu;,V2^^j^(^)^2D(^,^,t) into the GPE, Eq. (6a), after 
rescaling, we have 



where now 



-vi + V2 



2D 



(Bl) 



} 



/ dzwfj^{z) = 1/ ^/27ry. In the following we only present the 
integration of the last term of Eq. (Bl), 

dzdnn^{x, y, z, t)wfj^{z). (B3) 

) 

Again we use the identity — V^[/3D(r) = S{t) to write 

dzz {U3D * |V^2D^lDp) = -|V^2D^lDp 

- vi {U3D * |V^2D^lDp) • 

We recall that |V^2d(^, ^)^id(^)P is even in z so that 
U^^ * |V^2D^iDp also becomes even in z and d^^ziU^^ * 
|V^2D^iDp) becomes odd in z. Plugging ^nn = ^n^n^ + 
'Ti'ldzz + "^^zdn^z into Eq. (B3) yields 



3/^20 Cdd 



2D I 



dz {UsD * |V^2d^^^id|^) ^id(^) 



(B4) 



where the d^^z term disappears because the integral of odd 
functions vanishes. The first term in Eq. (B4) contributes to 
the modified contact interaction term in Eq. (16). 

After expanding the convolution in the last term of Eq. (B4), 
wee see that we need to compute the integral 



dz' dz 



Ati^{x - x'Y + - y'f + (z - z'Y 



(B5) 



By changing the variables according to z 
and integrating over z' we reduce Eq. (B5) to 



z — z', z' = z-\-z'. 



dz- 



1 



-52/27 



(B6) 



'2v^(2^)3/2 ^{x-x')^^{y-y')^^z^' 

By substituting p = z/^ and introducing polar coordinates 
with r = {x — x'Y^ {y — y'Y we get 



(27r)3/2 j_, 



dp 



-pV2 



Y^r^ + 7p2 



(B7) 



^ = -3/32Dedd j ^^V/73D(r-rO|(V^2D^lD)(r^^)|^ (B2) 

with rescaled ^id(^) = ( e"^, = d^^ + S^^^^, 

= /^2d(1 — Multiplying both sides of Eq. (Bl) by 

wiY){z) and integrating over the z direction, we can obtain a 
2D wave equation for ?/;2d. 

As in the preceding appendix, the integration is straightfor- 
ward for all but the last term because J dzw\^{z) = 1 and 



This integral representation of /72d(^) is identical to 
Eq. (17) [50, see Sec. 3.364]. We obtain the potential 
<I>2D, Eq. (16b), from the last term in Eq. (B4) after insert- 
ing Eq. (B7). Finally, plugging Eq. (B4) into the integrated 
Eq. (Bl) results in the mean-field equation for the quasi-2D 
dipolar BEC, Eq. (16). 

The Fourier transform U2D{kxj ky) of the radially symmet- 
ric kernel U2D is also radially symmetric and given by 



U2D{\kr 



1 

2^ 



dhi- 



\kr\ 



(B8) 



with \kr\ 



kx^k'^. The asymptotic behavior of this 



Fourier transform is t/2D(|^r|) 

UZ{\kr\) 



2T^\kr\ 

for \kr \ 00. 



for \kr 



and 
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Appendix C: Closed forms of the ID and 2D nonlocal potentials 



In this appendix, we derive closed forms for the nonlocal 
potentials <I>id, Eq. (9b), and $2D, Eq. (16b). Owing to the 
properties of the convolution in <I>id we may write Eq. (9b) as 



^iD = — ^ /-- (1 - 3n^ 



27r 



)(9..t/lD)*|V^lD|', 



i.e., the derivatives only affect the kernel Uij). Straightforward 
calculation of the second derivative of Eq. (10) leads to 



UM = d,,UM = 7^id(^)(1 + 7^') - V ^1^1- 



The resulting nonlocal potential $id is of the same form as the 
one given in Ref . [27] for the special case of axial polarization. 



In an analogous fashion, in Eq. (16b) is given by 



^2D = 



[(Sn^n^ -n^V')/72D] *|V^2d|'. (CI) 



By using Eq. (17) and the properties of the derivatives of 
Bessel functions we find for the second derivatives of U2D 



dyyU2D 



or' 1^1 



3" 

2(27r7)2 



+ — 

7 
7 

3^ 

2(27r7)2 



2x' 



2y^ 



xy 2xy\ 



Here and in the following we suppress the argument r^/47 of 



the Bessel functions with 
(CI) 



-\-y . Plugging these identi- 
ties into Eq. (CI) we find for the nonlocal kernel 1/2^^ (x, y) = 



2(27r7)3/2 



l-3n: 



{n^x + Uyyy 



{n^x + nyy)'^[l - 27/r^ 



(C2) 



and the nonlocal potential <I>2d = -^^^f^U^^^ * |V^2Dp. The 

kernel 1/2^^ (x, y) changes from a symmetric peak for = 1 
to a saddle shape with minima along the projection of the 
dipole axis onto the x-y plane for < 1. We note that 



^2D reduces to the potential for axial polarization (n^ = 1) 
derived, e.g., in Ref. [31]. However, Eq. (C2) is valid for arbi- 
trary polarization. 
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